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Abstract of EP0948237 

The noise reduction method uses subtraction of a 
noise compensation signal from the microphone 
signal (y), the noise compensation signal 
provided by simulation of the noise signal via an 
adaptive filter (H) using a least mean square 
algorithm, supplied with a reference signal, e.g. a 
loudspeaker signal, with transformation of the 
microphone signal, the compensation signal and 
the output signal into the frequency range. 



-.-™ ------ 



s. 



i 



Data supplied from the esp@cenet database - Worldwide 



http://v3.espacenet.com/textdoc?DB=EPODOC&IDX=EP0948237&F=0 



2007-07-3 1 



(19) 




Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 






(11) 



EP 0 948 237 A2 



(12) 



EUROPAISCHE PATENTANMELDUNG 



(43) VerOffentlichungstag: 

06.10.1999 Patentblatt 1999/40 

(21) Anmeldenummer: 99106123.5 

(22) Anmeldetag: 01.04.1999 



(51) Int. CI. 6 : H04R 3/00 



(84) Benannte Vertragsstaaten: 


(71) Anmelder: 


AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 


DaimlerChrysler Aerospace AG 


MCNLPTSE 


81663 Munchen (DE) 


Benannte Erstreckungsstaaten: 


(72) Erfinder: Thomas, Hans-Jorg 


AL LT LV MK RO SI 




89233 Neu-Ulm (DE) 


(30) Prioritat: 03.04.1998 DE 19814971 





CM 
< 

CO 
CM 

00 



(54) Verfahren zur Storbefreiung eines Mikrophonsignals 



(57) Fur ein Verfahren zur Befreiung eines Mikro- 
phonsignals von Storsignalanteilen durch Erzeugen 
eines Kompensationssignals und Subtraktion des Kom- 
pensationssignals vom Mikrophonsignal wird vorge- 
schlagen, die Kompensation vollstandig im 
Frequenzbereich vorzunehmen und auch das Aus- 
gangssignal im Frequenzbereich zu verarbeiten. Es 
werden MaBnahmen zur Verringerung des Aufwands 
bei der Signalverarbeitung angegeben. Vorteilhafte 
Weiterbildungen sehen beispielsweise vor, daB zur 



Storbefreiung eines Sprachsignals eine Filtereinstel- 
lung verwandt wird, welche in einer vorhergehenden 
Sprachpause gewonnen wurde, und/oder das fur lange 
Impulsantworten das Nachbildungsfilter in mehrere Teil- 
filter aufgespalten wird. Die Erfindung ist insbesondere 
vorteilhaft fur die Befreiung eines Spracheingabesi- 
gnals in einem Fahrzeug von Storsignalanteilen eines 
Radios Oder dgl., dessen Quellsignal als Referenzsi- 
gnal zur Verfugung steht. 
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Beschreibung 

[0001] Die Erfindung betrifft ein Verfahren zur Storbe- 
freiung eines Mikrophonsignals. 

[0002] Derartige Verfahren gewinnen insbesondere s 
fur die Spracheingabe von Kommandos und/oder fur 
Freisprechtelefone zunehmend an Bedeutung, wobei 
insbesondere die Situation in einem Fahrzeug ein wich- 
tiger Anwendungsfall ist. 

[0003] Eine besondere Situation ist in Fahrzeugen 10 
haufig dadurch gegeben, daB ein Wiedergabegerat wie 
z.B. ein Radio, ein Kassetten- oder CD-Abspieler uber 
einen Lautsprecher eine Gerauschumgebung erzeugt, 
die als Storsignal ein von einem Mikrophon aufgenom- 
menes Sprachsignal, beispielsweise fur die Spracher- is 
kennung oder Telefoniibertragung uberlagert. Zur 
Erkennung von Spracheingaben in einem Spracherken- 
ner oder zur verstand lichen SprachCibertragung uber 
Telefon ist das Mikrophonsignal soweit wie moglich von 
Storsignalanteilen zu befreien. 20 
[0004] Das von einer Storquelle, insbesondere einem 
Lautsprecher ausgehende Storsignal gelangt nicht nur 
auf direktem kurzestem Weg zum Mikrophon, sondern 
tritt auch noch uber zahlreiche Reflexionen als eine 
Uberlagerung einer Mehrzahl von Echos mit verschie- 25 
denen Laufzeiten im Mikrophonsignal in Erscheinung. 
Die gesamte Einwirkung des Storsignals von der Stor- 
quelle auf das Mikrophonsignal kann durch eine a priori 
unbekannte Ubertragungsfunktion des Raumes, bei- 
spielsweise des Fahrgastraumes eines Kraftfahrzeugs 30 
beschrieben werden. Die Ubertragungsfunktion andert 
sich je nach Besetzung des Fahrzeugs und nach Posi- 
tion der einzelnen Personen. Durch Nachbildung dieser 
Ubertragungsfunktion und Filterung eines Referenzsi- 
gnals von der Storquelle mit dieser Nachbildung kann 35 
ein Kompensationssignal erzeugt werden, welches 
durch Subtraktion vom Mikrophonsignal ein vom Storsi- 
gnal befreites Signal, beispielsweise ein reines Sprach- 
signal liefert. Im Realfall stellt die genannte 
Nachbildung eine mehr oder minder gute Annaherung 40 
an die unbekannte Ubertragungsfunktion dar und die 
Storung kann nicht vollstandig beseitigt werden. 
[0005] Aufgabe der vorliegenden Erfindung ist es, ein 
Verfahren zur Storbefreiung eines Mikrophonsignals 
anzugeben, daB bei vertretbarem Signalverarbeitungs- 45 
aufwand gute Eigenschaften hinsichltich der Entstorung 
aufweist. 

[0006] Die Erfindung ist im Patentanspruch 1 
beschrieben. Die Unteranspruche enthalten vorteilhafte 
Ausgestaltungen und Weiterbildungen der Erfindung. so 
[0007] Wesentlich an dem erfindungsgemaBen Ver- 
fahren ist, daB die Kompensation des Storsignalanteils 
im Mikrophonsignal mittels eines aus dem Referenzsi- 
gnal uber die Nachbildung der Ubertragungsfunktion 
erzeugten Kompensationssignals im Frequenzbereich 55 
vorgenommen wird, so daB Mikrophonsignal, Kompen- 
sationssignal und Ausgangssignal im Frequenzbereich, 
d.h. in Form von Spektren vorliegen. Die Signalverar- 



beitung in diesem Verfahrensschritt im Frequenzbereich 
erfordert zwar eine spektrale Transformation des Mikro- 
phonsignals, berucksichtigt aber, daB die Nachbildung 
der Ubertragungsfunktion im Frequenzbereich vorteil- 
hafter ist und stellt fur eine vorteilhafte nachfolgende 
zusatzliche Gerauschreduktion des Ausgangssignals, 
die typischerweise gleichfalls im Frequenzbereich vor- 
genommen wird, bereits eine besonders geeignete 
Signalform bereit. 

[0008] Durch einfache Naherungen beim Ersatz eines 
Verarbeitungsschritts mit einem Zeitfenster kann durch 
Ubergang zu einer Faltung im Frequenzbereich eine 
deutliche Reduzierung des Verarbeitungsaufwands 
reduziert werden. 

[0009] Fur lange Impulsantworten der Ubertragungs- 
funktion bzw. deren Nachbildung sieht eine vorteilhafte 
Weiterbildung der Erfindung eine Aufteilung des Nach- 
bildungsfilters in mehrere Teilfilter zu zeitversetzten 
Segmenten des segmentierten Referenzsignals vor, 
deren Koeffizienten-Aktualisierung zeitlich gestaffelt 
sein kann, wodurch der Signalverarbeitungsaufwand 
gering gehalten werden kann. 

[0010] Als besonders vorteilhaft erweist es sich, die 
Entstorung eines Sprachsignals auf der Basis einer Ein- 
stellung des Nachbildungsfilters, die in einer vorherge- 
henden Sprachpause gewonnen und gespeichert 
wurde, vorzunehmen. 

[001 1 ] Die Aufteilung des Nachbildungsfilters in meh- 
rere Teilfilter und die Storbefreiung auf der Basis einer 
in einer Sprachpause gewonnenen Filtereinstellung 
sind auch unabhangig von der Storsignalkompensation 
im Frequenzbereich eigenstandig fur die Storbefreiung 
eines Mikrophonsignals realisierbar und vorteilhaft. 
[0012] Die Erfindung ist nachfolgend anhand von 
bevorzugten Ausfuhrungsbeispielen unter Bezugnahme 
auf die Abbildungen noch eingehend veranschaulicht. 
Dabei zeigt: 

Fig. 1 ein Prinzip der Kompensation eines Radio- 
signals 

Fig. 2a ein Blockschaltbild zu Fig. 1 

Fig. 2b ein Blockschaltbild zur Filternachbildung 

Fig. 3 ein detailliertes Beispiel zu Fig. 2b 

Fig. 4 eine Erweiterung auf mehrere Teilfilter 

Fig. 5 einen Ubergang zur Kompensation im Fre- 
quenzbereich 

Fig. 6 ein detailliertes Beispiel zu Fig. 5b 

Fig. 7 ein Ausfiihrungsbeispiel mit mehreren Teil- 
filtern 

Fig. 8 ein Ausfiihrungsbeispiel mit Speicherung 
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der Filtereinstellungen 

Fig. 9 Signale einer synthetischen Beispielsszene 

Fig. 1 0 Impulsantwort und Ubertragungsf unktion zu 5 
Fig. 9 

Fig. 1 1 Signal einer ersten MeBszene 

Fig. 1 2 Impulsantwort und Ubertragungsf unktion zu 10 
Fig. 1 1 

Fig. 13 das Beispiel nach Fig. 11 mit Speicherung 
der Filtereinstellungen 

15 

Fig. 14 eine Sprachpausendetektion zu Fig. 13 

Fig. 15 Impulsantworten und Ubertragungsfunktio- 
nen zu Fig. 1 1 und Fig. 13 

20 

Fig. 16 Ubergang von einem Zeitfenster zu einer 
Faltung im Frequenzbereich 

Fig. 17 ein Rechteck-Zeitfenster mit Linienspek- 

trum 25 

Fig. 18 ein Hamming-Zeitfenster mit Linienspek- 
trum 

Fig. 19 Staffelung von Signalblocken bei der Filter- 30 
berechnung 

Fig. 20 Signale einer zweiten MeBszene 

Fig. 21 eine Sprachpausendetektion zu Fig. 20 35 

Fig. 22 Impulsantworten und Ubertragungsfunktio- 
nen zu Fig. 20 und Fig. 21 

Fig. 23 Signale einer dritten MeBszene 40 

Fig. 24 eine Sprachpausendetektion zu Fig. 23 

Fig. 25 Impulsantworten und Ubertragungsfunktio- 

nen zu Fig. 23 und Fig. 24 45 

Fig. 26 Signale einer vierten MeBszene 

Fig. 27 eine Sprachpausendetektion zu Fig. 26 

50 

Fig. 28 Impulsantworten und Ubertragungsfunktio- 
nen zu Fig. 26 und Fig. 27. 

[0013] Fig. 1 stellt das Prinzip einer Einrichtung zur 
(einkanaligen) Radiosignalkompensation dar. Das vom 55 
Lautsprecher abgestrahlte akustische Signal gelangt 
auf direktem Wege, aber auch uber zahlreiche Reflek- 
tionen im Fahrzeuginnenraum, auf das Mikrophon des 



Spracheingabesystems. Unter der Annahme, daB sich 
die Ubertragungsstrecke G demnach als Transversalfil- 
ter mit einer gewichteten Summe zeitlich verzogerter 
Echos darstellt, laBt sich eine Filternachbildung H fin- 
den, die im Idealfall H=G eine vollstandige Kompensa- 
tion des Radiosignales ermoglicht. 

[0014] Das Lautsprechersignal x wird durch die a 
priori unbekannte Ubertragungsfunktion G des Fahr- 
zeuginnenraumes gefiltert. Es entsteht die Storkompo- 
nente r, die sich mit dem Sprachsignal s zu dem 
Mikrophonsignal y addiert. Um die Storkomponente r zu 
kompensieren, wird mittels der Filternachbildung H ein 
Schatzwert r A aus dem Lautsprechersignal x erzeugt. 
Der Ausgang der Schaltung liefert den Schatzwert fur 
das Sprachsignal: 

s A = s + r - r A = s+E 

[0015] Dem Sprachsignal s ist also am Ausgang der 
Schaltung noch das Fehlersignal E=r-r A uberlagert, 
welches in der Praxis moglichst klein gehalten werden 
sollte. Das Sprachsignal kann noch Storungen in Form 
von z.B. Motorgerauschen oder externen Gerauschen 
enthalten, die aber in diesem Zusammenhang nicht 
explizit behandelt werden. 

[0016] H ist ein adaptives Filter und arbeitet nach 
einem in der Literatur bekannten Standard verfahr en, 
dem LMS-Algorithmus (least mean squares). Neben 
dem Eingangssignal x wird noch das Fehlersignal E 
benotigt, um die Koeffizientenadaption im Filter H zu 
bewerkstelligen. Hierfur ist das Ausgangssignal s A der 
Bestimmung der Filterkoeffizienten zugefuhrt. 
[0017] Fig. 2a zeigt in anderer Darstellung nochmals 
die Anordnung von Fig. 1 als Radiosignalkompensation. 
Das adaptive System H kann z.B. im Zeitbereich als 
FIR-Filter (finite-impulse-response-Filter) realisiert wer- 
den. Bei groBen Impulsantwortlangen, wie sie in der 
Praxis haufig auftreten, ist hierzu allerdings ein sehr 
hoher Rechenaufwand notwendig. Verschiedene Vor- 
teile gegenuber einer Zeitbereichsldsung bietet die 
Realisierung des LMS-Algorithmus im Frequenzbereich 
(FLMS). Wegen der blockweisen Verarbeitung von 
Daten in den als diskreten Fouriertransformationen rea- 
lisierten spektralen Transformationen und der Filterrea- 
lisierung im Frequenzbereich durch Multiplikationen 
wird dieses Verfahren besonders rechenzeitgunstig. 
[0018] Fig. 2b zeigt ein Blockschaltbild des FLMS- 
Algorithmus. Die zugehorige Theorie ist an sich 
bekannt und daher an dieser Stelle nicht im Detail 
behandelt. Es bedeuten F eine spektrale Transforma- 
tion FFT eines Zeitsignals in den Frequenzbereich und 
F" 1 die inverse IFFT. Die als Projektionen P1 , P2 und P3 
bezeichneten Verarbeitungsschritte dienen der korrek- 
ten Segmentierung der Daten durch die blockweise Ver- 
wendung mit der FFT bzw. IFFT und werden spater 
noch genauer erlautert. Die Arbeitsweise des Filters 
besteht in der Multiplikation des Referenzspektrums X 
mit dem Filter-Koeffizientenvektor H. Das Spektrum des 



3 



5 



EP 0 948 237 A2 



6 



Filterausgangs R A wird uber F" 1 zuriick in den Zeitbe- 
reich transformiert. Nach Anwendung der Projektion P2 
auf den Realteil des so erhaltenen Kompensationssi- 
gnals steht das Signal r A zur Verfugung. Die Differenz 
der Signale 5 

s A = y-r A = s + r-r A = s + E 

stellt den eigentlichen Ausgang, eine Schatzung der 
Spracheingabe, dar. 10 
[0019] Wesentlicher Bestandteil des adaptiven Filters 
ist die Koeffiiientenadaption im Block K, die im Fig. 2b 
durch die Erneuerungsgleichung 

H' = H' + AH' 15 

beschrieben wird. Die hier mit zwei spektralen Transfor- 
mationen besonders aufwendige Projektion P1 berech- 
net aus H' den fur die Filterung benotigten 
Koeffizientenvektor H. Zur Berechnung des Korrektur- 20 
vektors AH' wird neben dem Referenzspektrum X das 
Spektrum S A des mit P 3 bewerteten Ausgangssignales 
s + r - r A benotigt. 

[0020] Ein detailliertes Blockschaltbild des in Fig. 2b 
dargestellten FLMS-Algorithmus zeigt Fig. 3. Die Abt- 25 
astwerte eines Signals und die Stutzstellen der FFT 
seien in gebrauchlicher Weise als samples bezeichnet. 
Alle Spektraltransformationen und deren Inverse sind 
als 256-Punkte - FFT's, die jeweils um 128 samples 
iiberlappen, zu segmentieren. Zu beachten ist, daB sich 30 
das Ausgangssignal s A im Zeitbereich aus 128-sample- 
Blocken zusammensetzt. Es entsteht aus der Differenz 
der zweiten Blockhalften (also jeweils der samples 129 
bis 256) von Mikrofonsignal und gefiltertem Kompensa- 
tionssignal r A . Aufwendig ist die Projektion P1, die 2 35 
FFT's benotigt und den Vektor H' in den Vektor H 
umrechnet. Hierbei wird aus dem komplexen 256- 
Punkte-Ergebnisvektor der Rucktransformation vom 
Frequenz-in den Zeitbereich (IFFT) die erste Halfte 
(samples 1 bis 128) ausgeschnitten und die zweite 40 
Halfte (samples 129 bis 256) zu Null gesetzt. Nach 
Anwendung dieses Rechteckfensters im Zeitbereich 
erfolgt wieder mittels FFT die Transformation in den Fre- 
quenzbereich. Einfach ist die Projektion P2. Sie besteht 
aus der oben schon beschriebenen Ausschnittsbildung 45 
der letzten 128 samples, wodurch aus uberlappenden 
256-sample-Blocken wieder nicht uberlappende 128- 
sample-Blocke entstehen. Ebenfalls sehr einfach ist 
schlieBlich auch die Projektion P3, welche umgekehrt 
aus nicht uberlappenden 128-sample-Blocken des Aus- so 
gangssignales durch Voranstellen von 128 Nullwerten 
wiederum uberlappende 256-sample-Blocke bereit- 
stellt. Die Adaption der Filterkoeffizienten H' L+1 fur 
einen Zyklus L+1 besteht aus der Addition eines 
Erneuerungsvektors AH' L zum alten Koeffizientenvektor 55 
H' L . Diese Erneuerung errechnet sich aus dem Produkt 
zwischen dem Spektrum S A L des Ausgangssignals und 
dem konjugiert komplexen Spektrum X* L des Referenz- 



signales - gewichtet mit einer spektralen Leistungsnor- 
mierung 2ja L , AH' L =2|i L • X* L • S A L . Zum Zweck 
dieser Leistungsnormierung ist der mit einer Konstan- 
ten 2a multiplizierte Kehrwert des geglatteten Referenz- 
leistungsspektrums S xxL zu berechnen 
2ji L = 2a/S xx L , wozu ein rekursives Filter 1. Ordnung 
mit einer Konstanten p dient 

Sxx,L = P ' |X L | 2 + (1"P) ' S^. 

[0021] Die Arbeitsweise des LMS-Algorithmus wird 
erheblich von der Adaptionskonstante a und der Glat- 
tungskonstante p beeinfluBt. Zwischenspeicher in 
Rekursionsschleifen sind mit Sp bezeichnet. 
[0022] Die bisher beschriebene Anordnung des 
FLMS-Algorithmus erlaubt Filternachbildungen mit 
einer maximalen Impulsantwortlange von einer halben 
FFT-Lange, im Beispielsfall also 128 samples. Sollen 
langere Impulsantworten kompensiert werden, ist der 
schon bekannte FLMS-Algorithmus fur einen Teilf ilter 
(Fig. 4a) auf n Teilf ilter zu erweitern. Eine 3-Teilf ilter- 
Losung mit einer Impulsantwortlange von 3*128=384 
samples hat sich bei der Radiosignalunterdruckung im 
Pkw mit einem Spracheingabesystem bewahrt (Fig. 
4b). Der im Fig. 4a mit B bezeichnete Block mit den Ein- 
gangssignalen X und S A und dem Kompensations- 
Spektrum R A als Ausgang ist durch die im Fig. 4b dar- 
gestellte Erweiterung zu ersetzen. Das Spektrum X des 
Referenzsignals wird durch Zwischenspeicher D um 1 
bzw. 2 Blocklangen verzogert und das unverzogerte X1 
und die beiden verzogerten Spektren X2, X3 werden 
separat in mit in einer erweiterten Projektion P1 
getrennt bestimmten Koeffizientenvektoren H1, H2, H3 
multipliziert. Die Bildung der Koeffizientenvektoren 
erfolgt analog zum Fall nur eines Teilf ilters, wobei in K1 , 
K2, K3 jeweils das zugehorige Referenzspektrum mit 
dem Spektrum S A des Ausgangssignals verknupft wird. 
Der Aufwand wird hauptsachlich durch die Verdreifa- 
chung der Projektion P1 betrachtlich erhoht. Zusatzli- 
cher Speicherplatzbedarf wird notwendig um die 
Spektren des um 1 bzw. 2 Blocklangen zeitlich alteren 
Referenzsignales X bereitzustellen. 
[0023] Bei der beispielhaft angegebenen Aufgaben- 
stellung der Unterdruckung des Radiosignales bei 
Spracheingabe im Kfz ist es vorteilhaft die Ausgangsda- 
ten nicht im Zeit- sondern im Frequenzbereich auszuge- 
ben, da dadurch eine verbesserte Anpassung an eine 
nachgeschaltete Gerauschunterdruckung erreicht wer- 
den kann. Der bereits vorgestellte FLMS-Algorithmus 
mit einem Teilf ilter benotigt gemaB Fig. 5a insgesamt 5 
FFT's bei einem Ausgangssignal im Zeitbereich. Wird 
dem Ausgang eine FFT nachgeschaltet, erhoht sich der 
Aufwand bei einem Frequenzbereichs-Ausgangssignal 
auf 6 FFT's. Die gleiche FFT-Anzahl ergibt sich 
zunachst auch bei einer aquivalenten Losung nach Fig. 
5b. Diese Variante besitzt jedoch folgende Vorteile: 

Bei der zeitgleichen Spektralanalyse der Signale x 
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und y ist nur eine einzige 256-Punkte-FFT mit 
geringem Zusatzaufwand fur eine spektrale Sepa- 
ration notwendig. Man erzielt eine Einsparung von 
1 FFT. 

5 

Die hier mit P4 gekennzeichnete und neu def inierte 
Projektion ist bis auf das verwendete Zeitfenster 
formal identisch mit der Projektion P1. Wie spater 
gezeigt wird, laBt sich P4 durch eine relativ einfa- 
che Faltungsoperation im Frequenzbereich ersetz- 10 
ten, ohne daB eine merkliche EinbuBe an Qualitat 
in Kauf genommen werden muB. Man erzielt eine 
Einsparung von 2 FFT's. 

[0024] Fig. 6 stellt ein ausfuhrlicheres Blockschaltbild 15 
des FLMS-Algorithmus mit Frequenzbereichs-Aus- 
gangssignal dar und gestattet nochmals einen Ver- 
gleich mit Fig. 3 (Zeitbereichs-Ausgang). Unverandert 
geblieben ist die Filteradaption bestehend aus Glattung 
der spektralen Leistung, Leistungsnormierung und 20 
Koeffizientenerneuerung. Neu sind die FFT im Mikro- 
phonkanal, die Differenzbildung Y-R A im Frequenz- 
anstatt im Zeitbereich zur Ausgangsbildung, und 
schlieBlich die neu definierte Projektion P4, welche sich 
lediglich durch das komplementare Zeitbereichsfenster 25 
von der Projektion P1 unterscheidet. 
[0025] Als Vorstufe einer nachfolgend beschriebenen 
bevorzugten Ausfiihrung betrachte man Fig. 7. Darge- 
stellt ist der FLMS-Algorithmus mit 3 Teilfiltern (384- 
sample-lmpulsantwort), welcher eine ausreichende 30 
Unterdriickung des Radiosignales im Mikrophonkanal 
des Spracheingabesystems liefert. Die Projektionen P1 
und P4 sind vereinfacht dargestellt. Es ist der schon aus 
Fig. 4b bekannte zusatzliche Aufwand in Form der Spei- 
cher P sowie die Verdreifachung der Projektion P1 35 
ersichtlich. Im Gegensatz zur 1 -Teilf ilter-Losung nach 
Fig. 6 wird die Summe W von gegenwartigem und den 
beiden zeitlich vorausgegangenen Referenz-Leistungs- 
spektren auf den Eingang des rekursiven Filters gege- 
ben. Die Tatsache, daB am Filterausgang nun praktisch 40 
die 3-fache geglattete spektrale Leistung vorliegt, wird 
nach der Kehrwertbildung durch Multiplikation mit der 
Konstanten 6a berucksichtigt. Nach der spektralen Lei- 
stungsnormierung des in P4 modifizierten Ausgangs- 
spektrums S A wird die Filter-Adaption nun fur die 3 45 
Koeffizientenvektoren der 3 Teilf ilter separat durchge- 
fuhrt. 

[0026] Ein Beispiel Z0 fur die Wirkungsweise der 
Erfindung nach Fig. 7 zeigt Fig. 9. Die Eingangsdaten 
wurden synthetisch generiert. Das Referenzsignal X so 
stellt 100000 Abtastwerte eines weiBen GauBrau- 
schens bei einer Abtastfolgefrequenz von fs=12 kHz 
dar. Das Mikrophonsignal Y entstand durch Faltung die- 
ses Rauschsignales mit einer ebenfalls konstruierten 
384-sample-lmpulsantwort sowie der Addition eines 55 
auBerst schwachen Sprachsignales. Beim Abhoren die- 
ses in Fig. 9 oben aufgezeichneten Signals y sind die 10 
gesprochenen Ziffern gerade noch im farbigen (weil 



gefilterten) Rauschen zu erkennen. Das in den Zeitbe- 
reich zurucktransformierte Ausgangssignal des Schat- 
zers befreit nach einem ca. 1 Sekunde (12000 samples) 
dauernden Einschwingvorgang sehr wirkungsvoll die 
Spracheingabe vom Rauschen und liefert ein unverzerr- 
tes jedoch leicht verhalltes Sprachsignal S A (Fig. 9 
unten). Die verwendeten beiden Parameter waren 
<x=0,05 und p=0,5, Werte, die sich auch bei den spater 
vorgestellten Beispielen gut bewahrt haben. 

[0027] Aus den jeweils 1 29 samples langen Teil koeffi- 
zientenvektoren H1, H2, H3 der 3 Teilf ilter nach Bild7 
laBt sich nun zu jedem Zeitpunkt die hieraus resultie- 
rende 3*128-sample-lmpulsantwort bzw. die zugeho- 
rige Filterubertragungsfunktion berechnen. So zeigt 
Fig. 10 oben die 384-sample-lmpulsantwort, wie sie 
sich ganz am Ende der Szene- also nachdem die Ziffer 
"0" gesprochen wurde - ergibt. Sie ist ein sehr genaues 
Abbild derjenigen Impulsantwort, die zur Faltung mit 
weiBem GauBrauschen und damit zur synthetischen 
Generierung des Signales mikro benutzt wurde. Die 
zugehorige Betragsubertragungsfunktion (Fig. 10 
unten) im Bereich zwischen den Frequenzen 0 und fs/2 
= 6 kHz stellt einen mit zahlreichen schmalbandigen 
Resonanzuberhohungen behafteten TiefpaBfrequenz- 
gang dar. 

[0028] WeiBes Rauschen als Referenzeingangssignal 
und gefiltertes "farbiges" Rauschen als Mikrophonain- 
gangssignal stellen im Sinne der Aufgabenstellung, 
eine Nachbildung dieses Filters zu finden, den einfach- 
sten Fall dar. Da das Referenzsignal per Definitionen 
alle Frequenzanteile enthalt, gelingt die Filteradaption 
hier am schnellsten. Die zusatzliche additive Sprachein- 
gabe im Mikrophoneingangssignal - also das eigentli- 
che Nutzsignal des Spracheingabesystems - stellt fur 
den (F)LMS-Algorithmus eine Storung dar, welche die 
korrekte Adaption der Filterkoeffizienten behindert. 
Anders ausgedruckt: nur in Sprachpausen ist das 
System dazu in der Lage die Raumakustik des Fahr- 
zeuginnenraums (Strecke Radiolautsprecher zum 
Mikrophon) richtig nachzubilden und dadurch eine 
Kompensation der Radiowiedergabe zu bewirken. Im 
oben demonstrierten Beispiel gemaB Fig. 9 gelingt dies 
sehr gut, da der Mikrophon-eingang im wesentlichen 
aus Rauschen und nur zu einem sehr geringen Teil aus 
Spracheingabe besteht. 

[0029] Aus echten Messungen im Fahrzeug hingegen 
entstammten das an den Radiolautsprecherklemmen 
abgegriffene Referenzsignal radio und das vom Mikro- 
phon des Spracheingabesystem aufgezeichnete Signal 
mikro der Szene Z1. Dieses Mikrophonsignal ist in Fig. 
11 oben dargestellt, besteht aus 100000 samples und 
besitzt demnach bei einer Sampling-Frequenz von 12 
kHz eine zeitliche Dauer von ca. 8,3 Sekunden. Es han- 
delt sich um flussig und relativ schnelle gesprochene 
Sprache eines hinten rechts im Auto sitzenden Fahr- 
zeuginsassen, wahrend gleichzeitig Musik mit normaler 
Lautstarke aus dem Autoradio-Lautsprecher erklingt. 
Nach Anwendung der EntstormaBnahme gemaB Fig. 7 
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und Umrechnung in den Zeitbereich ergibt sich das in 
Fig. 1 1 unten dargestellte Ausgangssignal. Der Hortest 
ergibt eine deutliche Herausarbeitung des Sprachan- 
teils bzw. eine vor allem in den kurzen Sprachpausen 
bemerkenswerte Musikunterdriickung. Auffallig und von 
Nachteil ist jedoch, daB die erwunschte Radiosignalun- 
terdriickung in starkem MaBe davon abhangt, ob 
gerade gesprochen wird Oder nicht. Die wieder am Sze- 
nenende ermittelte 384-sample-lmpulsantwort mit 
zugehoriger Ubertragungsfunktion ist aus Fig. 12 
ersichtlich. Eine korrekte Impulsantwort ist an den typi- 
schen Nullsamples (Totzeit) am Anfang zu erkennen, 
welche von der Laufzeit des Direktschalls vom Radio- 
lautsprecher zum Mikrophon herruhren. Aus den hier 
vorhandenen starken Storungen am Anfang sowie am 
Ende der Impulsantwort laBt sich demnach der SchluB 
Ziehen, daB die Filteradaption an dieser Stelle wegen 
vorhandener Spracheingabe auBerst unzureichend ist. 

[0030] Die im folgenden anhand von Fig. 8 beschrie- 
bene Ausfuhrungsform beruht auf folgender Grundidee: 
ein geeignetes Merkmal dient zusammen mit einem 
Schwellenwert als Indikator fur eine Spracheingabe. 
Unterschreitet das Merkmal die Schwelle, so ist dies ein 
Anzeichen fur fehlende Spracheingabe. In diesem Fall 
kann - wie oben schon festgestellt - eine weitgehend 
ungestorte Filteradaption erfolgen. Bei Spracheingabe 
wird nun auf denjenigen Filterkoeffizientensatz zuruck- 
gegriffen, der unmittelbar vor der Schwellenuberschrei- 
tung - d.h. am Ende der vorangegangenen 
Sprachpause - abgespeichert wurde. Diese gespei- 
cherten Koeffizienten H10, H20, H30 liefern im Regelfall 
eine deutlich bessere Radiosignal-Kompensation als 
die unter dem storenden EinfluB der Spracheingabe 
sich standig andernden aktuellen Koeffizienten H, H2, 
H3. 

[0031] Fig. 8 stellt eine Ausfuhrung mit einer weiter 
verbesserten FLMS-Verarbeitung mit 3 Teilf iltern dar. 
Neben den schon in Fig. 7 vorhandenen aktuellen Fil- 
terkoeffizientenvektoren H1 , H2, H3, welche zur Bildung 
des fortlaufend adoptierten Ausgangssignals y-R beno- 
tigt wurden, existiert nun ein zusatzliches Ausgangssi- 
gnal (y-Ro), das unter Verwendung gespeicherter 
Koeffizienten H10, H20, H30 gebildet ist. Die aktuellen 
Koeffizientensatze H1, H2, H3 stellen nur bei fehlender 
Spracheingabe im eingeschwungenen Zustand ein 
brauchbares Kompensationsfilter im Frequenzbereich 
dar, liefern hingegen bei Spracheingabe ungenugende 
Filtereigenschaften, weil der AdaptionsprozeB in der 
Regelschleife standig gestort wird. Bei fehlender 
Spracheingabe d.h. hoher Filterqualitat sind die drei 
Schalter geschlossen und es werden die aktuellen 
Koeffizientensatze in die Koeffizientenspeicher M1, M2, 
M3 geschrieben: H10=H1, H20=H2, H30=H3. Die Aus- 
gange (y-Ro) und (y-Ra) sind identisch. Einsetzende 
Spracheingaben bewirken ein Offnen der 3 Schalter, 
wodurch die zuletzt in den Speichern M1, M2, M3 
befindlichen Koeffizienten H10, H20, H30 nicht mehr 
iiberschrieben werden und unverandert bleiben. Dieser 



Zustand, in welchem sich die Ausgange (Y-Ro) und (Y- 
Ra) unterschieden, wird solange beibehalten, bis wie- 
der eine Sprachpause detektiert und die Schalter 
geschlossen werden. 

5 [0032] Als Sprachpausenmerkmal fea hat sich die 
geglattete Summe aller Absolutwerte der Koeffizienten- 
korrekturvektoren AHV, AH2', AH3' bewahrt (Fig. 8a). 
Diese GroBe ist gleich Null bzw. weist kleine Zahlen- 
werte auf, wenn es keinen oder nur einen geringen 

10 Bedarf gibt, die Koeffizienten abzuandern. In Sprach- 
pausen ist dies der Fall, der Regelkreis ist praktisch ein- 
geschwungen. Storungen, wie sie durch 
Spracheingabe - aber auch durch Bewegungen der 
Fahrzeuginsassen - hervorgerufen werden, haben 

15 einen erhohten Nachregelbedarf zur Folge, was sich 
durch entsprechend groBe Zahlenwerte bei AHV, AH2', 
AH3' und somit beim Merkmal fea bemerkbar macht. 
Ein Glattungsf ilter beispielsweise ein rekursiver TiefpaB 
1 . Ordnung mit dem Eingang feat stellt an seinem Aus- 

20 gang das geglattete Sprachpausen-Merkmal fea zur 
Verfugung, welches nach Vergleich mit einem Schwell- 
wert th die Schalter fur die Koeffizientenubernahme 
steuert. 

[0033] Die Wirkungsweise des verbesserten FLMS- 

25 Algorithmus nach Fig. 8 demonstriert Fig. 13. Oben ist 
das aufgezeichnete Signal y der Szene Z1 (vgl. Fig. 1 1 
oben) dargestellt, unten das gewonnene Ausgangssi- 
gnal. Schon der visuelle Vergleich der Ausgangssignale 
von Fig. 13 und Fig. 11 zeigt die verbesserte Herausar- 

30 beitung der Sprachpassagen. Der vergleichende Hor- 
test bestatigt dies: auch wahrend der Spracheingabe ist 
die Musikunterdruckung deutlich besser. Den Verlauf 
des Sprachpausenmerkmals und der konstanten 
Schwelle uber der Zeit (hier in FFT-Blocken skaliert) 

35 zeigt Fig. 14 oben. In den durch die Schwellenunter- 
schreitung detektierten Sprachpausen (Fig. 14 unten) 
findet laufend die Ubernahme der Koeffizienten in die 
Speicher wie beschrieben statt, um dort wahrend der 
Spracheingabe als gespeicherte Koeffizienten zur Ver- 

40 fugung zu stehen. Die schon in Fig. 12 am Szenenende 
gemessene 384-sample-lmpulsantwort mit zugehoriger 
Betragsubertragungsfunktion ist in Fig. 15 als aktuelle 
Impulsantwort (a) bzw. aktuelle Ubertragungsfunktion 
(b) dargestellt. Im Gegensatz zu dieser infolge Sprach- 

45 eingabe stark gestorten Schatzung aus den aktuellen 
Koeffizienten H1, H2, H3 ist aus den gespeicherten 
Koeffizienten H10, H20, H30 eine Impulsantwort (c) und 
eine Ubertragungsfunktion (d) hoher Qualitat berechen- 
bar. Die Impulsantwort aus den gespeicherten Koeffizi- 

50 enten weist die typischen Nullsamples am Anfang auf, 
welche durch die Laufzeit des Direktschalls vom Radio- 
lautsprecher zum Spracheingabemikrophon verursacht 
werden. Aus der im Beispielsfall abzulesenden Totzeit 
von ca. 40 samples laBt sich die Entfernung zwischen 

55 Lautsprecher und Mikrophon bestimmen. 

[0034] Wie vorstehend schon angedeutet laBt sich die 
aufwendige Projektion P4 (IFFT, Fenster rechts im Zeit- 
bereich, FFT) ohne merkliche EinbuBe an Qualitat 
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durch eine relativ einfache Faltung im Frequenzbereich 
ersetzten, wodurch 2 FFT's eingespart werden. Dazu 
betrachte man Fig. 16. In einem ersten Schritt wird das 
"rechtsseitige" 128-sample-Rechteckfenster im Zeitbe- 
reich (Fig. 16a) bei der idealen Projektion ersetzt durch 5 
ein 128-sample-Hammingfenster (Fig. 16b). Gegenuber 
dem Rechteckfenster besitzt dieses den Vorteil eines 
bedeutend schmaleren Spektrums. Wie Fig. 17 zeigt, 
besteht beim Rechteckfenster der Realteil des Spek- 
trums aus einer einzigen Linie (Gleichanteil), wahrend w 
das zur Mitte antisymmetrische Imaginarteil-Spektrum 
aus vielen nach auBen hin langsam abfallenden Linien 
mit alternierenden Nullstellen besteht. Im Gegensatz 
dazu beschrankt sich das komplexe Spektrum des 
Hammingfensters (Fig. 18) auf insgesamt 7 Linien, von 15 
denen im symmetrischen Realteil nur 3 und im antisym- 
metrischen Imaginarteil nur 4 Werte von Null verschie- 
den sind. Samtliche weiter auBen liegenden Anteile 
sind vernachlassigbar gering. Diese spezielle Eigen- 
schaft des Hammingfensters ermoglicht es vorteilhaf- 20 
terweise die Multiplikation im Zeitbereich (Fig. 16b) zu 
ersetzen durch eine Faltung mit dem zugehorigen 7- 
sample-Spektrum im Frequenzbereich und damit eine 
IFFT und eine FFT einzusparen (Fig. 16c). 

[0035] Prinzipiell laBt sich naturlich auch die Projek- 25 
tion P1 (IFFT - linksseitiges Rechteckfenster - FFT) 
ersetzen durch eine entsprechende Faltungsoperation 
im Frequenzbereich mit dem konjugiert komplexen 7- 
Linien-Spektrum. Experimente haben jedoch gezeigt, 
daB Einsparungen an dieser Stelle erkauft werden mit 30 
einer deutlichen Verschlechterung des Einschwingver- 
haltens. Aufwandsgunstige Losungen lassen sich trotz- 
dem dadurch erzielen, daB in dem LMS-Algorithmus 
nach Fig. 8 die 3 Projektionen P1 nicht gleichzeitig in 
einem 256-sample-lnputdatenblock abgearbeitet wer- 35 
den mussen. Die mit 128-samples uberlappenden 
Inputdatenblocke der Lange 256 sind mit einer willkur- 
lich bei "1" beginnenden Nummerierung in Fig. 19a 
skizziert. So ist es z.B. moglich bei modulo-3-Zahlweise 
der Inputdatenblocke die 3 Teilfilterprojektionen nicht 40 
parallel (Fig. 19b) sondern sequentiell in aufeinanderfol- 
genden Blocken Fig. 19 zu berechnen. Dadurch sind 
bei idealer Projektion P1 pro Datenblock nicht 6 son- 
dern nur noch 2 FFT's notwendig. Es hat sich gezeigt, 
daB die Kompensation des Radiosignales auch noch 45 
ausreichend funktioniert, wenn die Abstande zwischen 
den zu berechnenden Teilfilter-Projektionen noch gro- 
Ber gewahlt werden. Zahlt man die Blocke z.B. modulo 
6, so ist lediglich in jedem zweiten Block eine Projektion 
Zu berechnen (Fig. 19d). Selbst eine Reduzierung auf so 
einen Abstand von vier Blocken zwischen zwei aufein- 
anderfolgenden P1 Berechnungen mittels modulo-12- 
Zahlung fuhrt noch zu brauchbaren Ergebnissen (Fig. 
19e). 

[0036] Die Leistungsfahigkeit des FLMS-Algorithmus 55 
mit 3 Teilf iltern gemaB Blockschaltung Fig. 8 und einer 
sequentiellen Berechnung der idealen Projektion P1 im 
Zeitraster nach Fig. 19e sowie der Projektion P2 mittels 



Faltung im Frequenzbereich (Fig. 16c) mit einem kom- 
plexen 7-Linien-Spektrum (Fig. 18) sei anhand von 3 
MeBszenen demonstriert. 

[0037] Die erste dieser Szenen Z2 beinhaltet Sprach- 
eingabe von Ziffern, wobei der Radiolautsprecher anna- 
hernd weiBes Rauschen mit verhaltnismaBig hoher 
Lautstarke abstrahlt. Das zugehorige 100000-sample- 
Mikrophonsignal ist in Fig. 20 oben, das extrahierte 
Ausgangssignal in Fig. 20 unten dargestellt. Eine deut- 
liche Rauschbefreiung des Outputsignales gegenuber 
dem Mi krophon input stellt man durch Abhorvergleich 
fest. Der zeitliche Verlauf des Sprachpausenmerkmals 
ist zusammen mit der konstanten Schwelle th Fig. 21 
oben abgebildet und die hieraus abgeleiteten Sprach- 
pausen bzw. die zugeordneten Schalterstellungen in 
Fig. 21 unten. SchlieBlich zeigt Fig. 22 in zu Fig. 15 ana- 
loger Weise die am Szenenende gefundene Impulsant- 
wort (a) und Ubertragungsfunktion (b) auf der Basis der 
aktuellen Koeffizienten und die entsprechenden Gro- 
Ben (c), (d) auf der Basis der Sprachpauseneinstellung. 
Es ist deutlich erkennbar, daB die am Szenenende 
gefundene aktuelle Impulsantwort ein infolge Sprach- 
eingabe gestortes Ergebnis darstellt, wahrend die aus 
der letzten Sprachpause stammende Impulsantwort 
aus den gespeicherten Koeffizientensatzen eine hohe 
Qualitat aufweist. 

[0038] Die ersten 100000 samples einer MeBszene 
Z3 mit POP-Musik im Radio und flussig bis schnell 
gesprochener Sprache der rechts hinten sitzenden Per- 
son sind in Form des Mikrophonsignales y in Fig. 23 
oben aufgezeichnet. Nach ca. 10000 samples (0,83 s) 
wird das Radiosignal brauchbar unterdriickt (Fig. 23 
unten). Auch bei der im letzten Drittel dieser Szene ein- 
setzenden Spracheingabe bleibt die POP-Musikunter- 
druckung wirksam erhalten, wodurch die 
Sprachverstandlichkeit hier gegenuber dem Mikrophon- 
signal merklich verbessert wird. Nach einer langen 
Sprachpause kommt es wegen der anschlieBenden 
pausenfreien Spracheingabe nicht mehr zu einer 
Schwellenunterschreitung (Fig. 24). Aus diesem 
Grunde ist die in Fig. 25 unten am Ende der Szene fest- 
gehaltene Impulsantwort auf der Basis der gespeicher- 
ten Koeffizienten zeitlich relativ veraltet, weil sie bereits 
ca. 2,3 svorher aktuell war (215 Blocke * 10,7 ms). Wie- 
der weist die aktuelle Impulsantwort (Fig. 25 oben) 
starke von der Spracheingabe herruhrende Storungen 
auf. Wie ein Vergleich mit der ahnlichen Szene Z1 nach 
Figuren 11 bis 15 zeigt, ist trotz des stark verringerten 
Rechenaufwandes die Qualitat der Storbefreiung unver- 
andert hoch. 

[0039] Die letzte Szene Z4 nach Fig. 26 wurde ohne 
Spracheingabe erstellt und soli abschlieBend nochmals 
die Musikunterdruckungseigenschaften des beschrie- 
benen FLMS-Algorithmus demonstrieren. Nach ca. 
18000 samples bzw. 1 ,5 s wird - wie aus Fig. 26 unten 
ersichtlich - die Musik wirksam unterdriickt. Diese 
Eigenschaft wird bis zum Szenenende mit unverander- 
ter Qualitat beibehalten. Fig. 27 zeigt auf, daB das 
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Sprachpausen-GroBe fea uberwiegend unter der 
Schwelle th bleibt. Die Zeiten, in welchen auf die 
gespeicherten Koeffizienten zuruckgegriffen wird, sind 
demnach nur sehr kurz. Impulsantwort und Ubertra- 
gungsfunktion aus aktuellen Koeffizienten sind daher im 5 
wesentlichen mit den entsprechenden Verlaufen aus 
Sprachpausen-Koeffizienten identisch. 

Patentanspruche 

10 

1. Verfahren zur Storbefreiung eines Mikrophonsi- 
gnals von Anteilen eines Quellsignals, daB als 
Referenzsignal (x) vorliegt und nach Durchlaufen 
einer Ubertragungsstrecke mit a priori unbekannter 
Ubertragungsfunktion (G) sich im Mikrophonsignal 15 
als Storsignal (r) einem Sprachsignal (s) uberla- 
gert, durch adaptive Nachbildung des Storsignals 
und Kompensation des tatsachlichen und des 
nachgebildeten Storsignals in einem Ausgangssi- 
gnal, wobei das Mikrophonsignal gleichfalls in den 20 
Frequenzbereich transformiert, die Kompensation 

an Signalen im Frequenzbereich vorgenommen 
und das im Frequenzbereich vorliegende Aus- 
gangssignal zur Adaption der Nachbildung mit dem 
im Frequenzbereich vorliegenden Referenzsignal 25 
verknupft wird. 

2. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet, daB das Ausgangssignalspektrum in den 
Zeitbereich transformiert, das Zeitsignal durch Vor- 30 
anstellen von Nullwerten auf doppelte Lange 
gebracht, in den Frequenzbereich rucktransformiert 
und der Nachbildung der Ubertragungsfunktion 
zugrundegelegt wird. 

35 

3. Verfahren nach Anspruch 1, dadurch gekenn- 
zeichnet, daB das Ausgangssignalspektrum mit 
dem Spektrum eines Hamming-Zeitfensters gefal- 
tet wird und der Nachbildung der Ubertragungs- 
funktion zugrundegelegt wird. 40 

4. Verfahren nach einem der Anspruche 1 bis 3, 
dadurch gekennzeichnet, daB zur Nachbildung des 
Storsignalanteils eine adaptive Filterfunktion eines 
Nachbildungsfilters auf das Referenzsignal ange- 45 
wandt wird. 

5. Verfahren nach Anspruch 4, dadurch gekenn- 
zeichnet, daB die Filterfunktion durch einen Koeffi- 
zientenvektor vorgegeben wird, dessen so 
Koeffizienten adaptiv eingestellt werden. 

6. Verfahren nach Anspruch 4 oder 5, dadurch 
gekennzeichnet, daB das Auftreten eines Sprachsi- 
gnalanteils im Mikrophonsignal detektiert wird und 55 
bei Auftreten eines Sprachsignals die vor Auftreten 
des Sprachsignals eingestellte Filterfunktion zur 
Bildung des Ausgangssignals beibehalten wird. 



7. Verfahren nach Anspruch 6, dadurch gekenn- 
zeichnet, daB auch bei Detektion eines Sprachsi- 
gnals die adaptive Nachfuhrung einer aktuellen 
Filterfunktion zusatzlich zur Bildung des Ausgangs- 
signals fortgefuhrt wird. 

8. Verfahren nach Anspruch 7, dadurch gekenn- 
zeichnet, daB das Auftreten eines Sprachsignals 
aus einer Veranderung der aktuellen Filterfunktion 
detektiert wird. 

9. Verfahren nach Anspruch 8, dadurch gekenn- 
zeichnet, daB die Veranderung der aktuellen Filter- 
funktion fur die Detektion des Auftretens eines 
Sprachsignals zeitlich geglattet wird. 

10. Verfahren nach einem der Anspruche 4 bis 9, 
dadurch gekennzeichnet, daB die Filterfunktion in 
mehrere Teilfilterfunktion zu aufeinanderfolgenden 
Abschnitten einer Gesamt- Impulsantwort aller Teil- 
filter aufgespalten ist und auf Referenzsignalspek- 
tren zu zeitlich versetzten Zeitsegmenten des 
segmentierten Referenz-Zeitsignals angewandt 
wird. 

11. Verfahren nach Anspruch 10, dadurch gekenn- 
zeichnet, daB die Adaption der Filterfunktion fur die 
Teilf ilter parallel durchgefiihrt wird. 

12. Verfahren nach Anspruch 10, dadurch gekenn- 
zeichnet, daB die Adaption der Filterfunktion fur die 
einzelnen Teilf ilter zeitsequentiell durchgefiihrt 
wird. 
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